Isolation of single circulating trophoblasts from maternal circulation for noninvasive fetal copy number variant profiling.

OBJECTIVE: To develop a multi-step workflow for the isolation of circulating extravillous trophoblasts (cEVTs) by describing the key steps enabling a semi-automated process, including a proprietary algorithm for fetal cell origin genetic confirmation and copy number variant (CNV) detection. METHODS: Determination of the limit of detection (LoD) for submicroscopic CNV was performed by serial experiments with genomic DNA and single cells from Coriell cell line biobank with known imbalances of different sizes. A pregnancy population of 372 women was prospectively enrolled and blindly analyzed to evaluate the current workflow. RESULTS: An LoD of 800 Kb was demonstrated with Coriell cell lines. This level of resolution was confirmed in the clinical cohort with the identification of a pathogenic CNV of 800 Kb, also detected by chromosomal microarray. The mean number of recovered cEVTs was 3.5 cells per sample with a significant reverse linear trend between gestational age and cEVT recovery rate and number of recovered cEVTs. In twin pregnanices, evaluation of zygosity, fetal sex and copy number profiling was performed in each individual cell. CONCLUSION: Our semi-automated methodology for the isolation and single-cell analysis of cEVTS supports the feasibility of a cell-based noninvasive prenatal test for fetal genomic profiling.

Circulating tumor cell copy-number heterogeneity in ALK-rearranged non-small-cell lung cancer resistant to ALK inhibitors.

Gatekeeper mutations are identified in only 50% of the cases at resistance to Anaplastic Lymphoma Kinase (ALK)-tyrosine kinase inhibitors (TKIs). Circulating tumor cells (CTCs) are relevant tools to identify additional resistance mechanisms and can be sequenced at the single-cell level. Here, we provide in-depth investigation of copy number alteration (CNA) heterogeneity in phenotypically characterized CTCs at resistance to ALK-TKIs in ALK-positive non-small cell lung cancer. Single CTC isolation and phenotyping were performed by DEPArray or fluorescence-activated cell sorting following enrichment and immunofluorescence staining (ALK/cytokeratins/CD45/Hoechst). CNA heterogeneity was evaluated in six ALK-rearranged patients harboring ≥ 10 CTCs/20 mL blood at resistance to 1st and 3rd ALK-TKIs and one presented gatekeeper mutations. Out of 82 CTCs isolated by FACS, 30 (37%) were ALK+/cytokeratins-, 46 (56%) ALK-/cytokeratins+ and 4 (5%) ALK+/cytokeratins+. Sequencing of 43 CTCs showed highly altered CNA profiles and high levels of chromosomal instability (CIN). Half of CTCs displayed a ploidy >2n and 32% experienced whole-genome doubling. Hierarchical clustering showed significant intra-patient and wide inter-patient CTC diversity. Classification of 121 oncogenic drivers revealed the predominant activation of cell cycle and DNA repair pathways and of RTK/RAS and PI3K to a lower frequency. CTCs showed wide CNA heterogeneity and elevated CIN at resistance to ALK-TKIs. The emergence of epithelial ALK-negative CTCs may drive resistance through activation of bypass signaling pathways, while ALK-rearranged CTCs showed epithelial-to-mesenchymal transition characteristics potentially contributing to ALK-TKI resistance. Comprehensive analysis of CTCs could be of great help to clinicians for precision medicine and resistance to ALK-targeted therapies.

Somatic Focal Copy Number Gains of Noncoding Regions of Receptor Tyrosine Kinase Genes in Treatment-Resistant Epilepsy.

Epilepsy is a heterogenous group of disorders defined by recurrent seizure activity due to abnormal synchronized activity of neurons. A growing number of epilepsy cases are believed to be caused by genetic factors and copy number variants (CNV) contribute to up to 5% of epilepsy cases. However, CNVs in epilepsy are usually large deletions or duplications involving multiple neurodevelopmental genes. In patients who underwent seizure focus resection for treatment-resistant epilepsy, whole genome DNA methylation profiling identified 3 main clusters of which one showed strong association with receptor tyrosine kinase (RTK) genes. We identified focal copy number gains involving epidermal growth factor receptor (EGFR) and PDGFRA loci. The dysplastic neurons of cases with amplifications showed marked overexpression of EGFR and PDGFRA, while glial and endothelial cells were negative. Targeted sequencing of regulatory regions and DNA methylation analysis revealed that only enhancer regions of EGFR and gene promoter of PDGFRA were amplified, while coding regions did not show copy number abnormalities or somatic mutations. Somatic focal copy number gains of noncoding regulatory represent a previously unrecognized genetic driver in epilepsy and a mechanism of abnormal activation of RTK genes. Upregulated RTKs provide a potential avenue for therapy in seizure disorders.

Genomic profiles of primary and metastatic esophageal adenocarcinoma identified via digital sorting of pure cell populations: results from a case report.

BACKGROUND: We report on a female patient who underwent primary radical resection for a stage 2B Her-2-positive Barrett's-type esophageal adenocarcinoma (EAC). Despite Her-2 targeted therapy, her disease recurred and required repeated metastectomies. CASE PRESENTATION: Digital cell sorting and targeted sequencing of cancer sub-clones from EAC and metastases revealed a completely mutated TP53, whereas the sorted stromal cells were wild-type. Her-2 amplification was significantly lower in the metastases when the patient became therapy-resistant. CONCLUSIONS: The mechanism of therapy resistance illustrated by this case could only be detected through accurate analysis of tumor sub-populations. Investigating tumor sub-populations of recurrent disease is important for adjusting therapy in recurrent EAC.

High-Throughput Sequencing of microRNAs in Glucocorticoid Sensitive Paediatric Inflammatory Bowel Disease Patients.

The aim of this research was the identification of novel pharmacogenomic biomarkers for better understanding the complex gene regulation mechanisms underpinning glucocorticoid (GC) action in paediatric inflammatory bowel disease (IBD). This goal was achieved by evaluating high-throughput microRNA (miRNA) profiles during GC treatment, integrated with the assessment of expression changes in GC receptor (GR) heterocomplex genes. Furthermore, we tested the hypothesis that differentially expressed miRNAs could be directly regulated by GCs through investigating the presence of GC responsive elements (GREs) in their gene promoters. Ten IBD paediatric patients responding to GCs were enrolled. Peripheral blood was obtained at diagnosis (T0) and after four weeks of steroid treatment (T4). MicroRNA profiles were analyzed using next generation sequencing, and selected significantly differentially expressed miRNAs were validated by quantitative reverse transcription-polymerase chain reaction. In detail, 18 miRNAs were differentially expressed from T0 to T4, 16 of which were upregulated and 2 of which were downregulated. Out of these, three miRNAs (miR-144, miR-142, and miR-96) could putatively recognize the 3’UTR of the GR gene and three miRNAs (miR-363, miR-96, miR-142) contained GREs sequences, thereby potentially enabling direct regulation by the GR. In conclusion, we identified miRNAs differently expressed during GC treatment and miRNAs which could be directly regulated by GCs in blood cells of young IBD patients. These results could represent a first step towards their translation as pharmacogenomic biomarkers.

A streamlined workflow for single-cells genome-wide copy-number profiling by low-pass sequencing of LM-PCR whole-genome amplification products.

Chromosomal instability and associated chromosomal aberrations are hallmarks of cancer and play a critical role in disease progression and development of resistance to drugs. Single-cell genome analysis has gained interest in latest years as a source of biomarkers for targeted-therapy selection and drug resistance, and several methods have been developed to amplify the genomic DNA and to produce libraries suitable for Whole Genome Sequencing (WGS). However, most protocols require several enzymatic and cleanup steps, thus increasing the complexity and length of protocols, while robustness and speed are key factors for clinical applications. To tackle this issue, we developed a single-tube, single-step, streamlined protocol, exploiting ligation mediated PCR (LM-PCR) Whole Genome Amplification (WGA) method, for low-pass genome sequencing with the Ion Torrent™ platform and copy number alterations (CNAs) calling from single cells. The method was evaluated on single cells isolated from 6 aberrant cell lines of the NCI-H series. In addition, to demonstrate the feasibility of the workflow on clinical samples, we analyzed single circulating tumor cells (CTCs) and white blood cells (WBCs) isolated from the blood of patients affected by prostate cancer or lung adenocarcinoma. The results obtained show that the developed workflow generates data accurately representing whole genome absolute copy number profiles of single cell and allows alterations calling at resolutions down to 100 Kbp with as few as 200,000 reads. The presented data demonstrate the feasibility of the Ampli1™ WGA-based low-pass workflow for detection of CNAs in single tumor cells which would be of particular interest for genome-driven targeted therapy selection and for monitoring of disease progression.

QueryOR: a comprehensive web platform for genetic variant analysis and prioritization.

BACKGROUND: Whole genome and exome sequencing are contributing to the extraordinary progress in the study of human genetic variants. In this fast developing field, appropriate and easily accessible tools are required to facilitate data analysis. RESULTS: Here we describe QueryOR, a web platform suitable for searching among known candidate genes as well as for finding novel gene-disease associations. QueryOR combines several innovative features that make it comprehensive, flexible and easy to use. Instead of being designed on specific datasets, it works on a general XML schema specifying formats and criteria of each data source. Thanks to this flexibility, new criteria can be easily added for future expansion. Currently, up to 70 user-selectable criteria are available, including a wide range of gene and variant features. Moreover, rather than progressively discarding variants taking one criterion at a time, the prioritization is achieved by a global positive selection process that considers all transcript isoforms, thus producing reliable results. QueryOR is easy to use and its intuitive interface allows to handle different kinds of inheritance as well as features related to sharing variants in different patients. QueryOR is suitable for investigating single patients, families or cohorts. CONCLUSIONS: QueryOR is a comprehensive and flexible web platform eligible for an easy user-driven variant prioritization. It is freely available for academic institutions at http://queryor.cribi.unipd.it/ .

Digital Sorting of Pure Cell Populations Enables Unambiguous Genetic Analysis of Heterogeneous Formalin-Fixed Paraffin-Embedded Tumors by Next Generation Sequencing.

Precision medicine in oncology requires an accurate characterization of a tumor molecular profile for patient stratification. Though targeted deep sequencing is an effective tool to detect the presence of somatic sequence variants, a significant number of patient specimens do not meet the requirements needed for routine clinical application. Analysis is hindered by contamination of normal cells and inherent tumor heterogeneity, compounded with challenges of dealing with minute amounts of tissue and DNA damages common in formalin-fixed paraffin-embedded (FFPE) specimens. Here we present an innovative workflow using DEPArray™ system, a microchip-based digital sorter to achieve 100%-pure, homogenous subpopulations of cells from FFPE samples. Cells are distinguished by fluorescently labeled antibodies and DNA content. The ability to address tumor heterogeneity enables unambiguous determination of true-positive sequence variants, loss-of-heterozygosity as well as copy number variants. The proposed strategy overcomes the inherent trade-offs made between sensitivity and specificity in detecting genetic variants from a mixed population, thus rescuing for analysis even the smaller clinical samples with low tumor cellularity.

A deep survey of alternative splicing in grape reveals changes in the splicing machinery related to tissue, stress condition and genotype.

BACKGROUND: Alternative splicing (AS) significantly enhances transcriptome complexity. It is differentially regulated in a wide variety of cell types and plays a role in several cellular processes. Here we describe a detailed survey of alternative splicing in grape based on 124 SOLiD RNAseq analyses from different tissues, stress conditions and genotypes. RESULTS: We used the RNAseq data to update the existing grape gene prediction with 2,258 new coding genes and 3,336 putative long non-coding RNAs. Several gene structures have been improved and alternative splicing was described for about 30% of the genes. A link between AS and miRNAs was shown in 139 genes where we found that AS affects the miRNA target site. A quantitative analysis of the isoforms indicated that most of the spliced genes have one major isoform and tend to simultaneously co-express a low number of isoforms, typically two, with intron retention being the most frequent alternative splicing event. CONCLUSIONS: As described in Arabidopsis, also grape displays a marked AS tissue-specificity, while stress conditions produce splicing changes to a minor extent. Surprisingly, some distinctive splicing features were also observed between genotypes. This was further supported by the observation that the panel of Serine/Arginine-rich splicing factors show a few, but very marked differences between genotypes. The finding that a part the splicing machinery can change in closely related organisms can lead to some interesting hypotheses for evolutionary adaptation, that could be particularly relevant in the response to sudden and strong selective pressures.

Chromosome scale genome assembly and transcriptome profiling of Nannochloropsis gaditana in nitrogen depletion.

Nannochloropsis is rapidly emerging as a model organism for the study of biofuel production in microalgae. Here, we report a high-quality genomic assembly of Nannochloropsis gaditana, consisting of large contigs, up to 500 kbp long, and scaffolds that in most cases span the entire length of the chromosomes. We identified 10646 complete genes and characterized possible alternative transcripts. The annotation of the predicted genes and the analysis of cellular processes revealed traits relevant for the genetic improvement of this organism such as genes involved in DNA recombination, RNA silencing, and cell wall synthesis. We also analyzed the modification of the transcriptional profile in nitrogen deficiency-a condition known to stimulate lipid accumulation. While the content of lipids increased, we did not detect major changes in expression of the genes involved in their biosynthesis. At the same time, we observed a very significant down-regulation of mitochondrial gene expression, suggesting that part of the Acetyl-CoA and NAD(P)H, normally oxidized through the mitochondrial respiration, would be made available for fatty acids synthesis, increasing the flux through the lipid biosynthetic pathway. Finally, we released an information resource of the genomic data of N. gaditana, available online at www.nannochloropsis.org.

PASS-bis: a bisulfite aligner suitable for whole methylome analysis of Illumina and SOLiD reads.

SUMMARY: The sequencing of bisulfite-treated DNA (Bi-Seq) is becoming a gold standard for methylation studies. The mapping of Bi-Seq reads is complex and requires special alignment algorithms. This problem is particularly relevant for SOLiD color space, where the bisulfite conversion C/T changes two adjacent colors into 16 possible combinations. Here, we present an algorithm that efficiently aligns Bi-Seq reads obtained either from SOLiD or Illumina. An accompanying methylation-caller program creates a genomic view of methylated and unmethylated Cs on both DNA strands. AVAILABILITY AND IMPLEMENTATION: The algorithm has been implemented as an option of the program PASS, freely available at http://pass.cribi.unipd.it. CONTACT: pass@cribi.unipd.it SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

First survey of the wheat chromosome 5A composition through a next generation sequencing approach.

Wheat is one of the world's most important crops and is characterized by a large polyploid genome. One way to reduce genome complexity is to isolate single chromosomes using flow cytometry. Low coverage DNA sequencing can provide a snapshot of individual chromosomes, allowing a fast characterization of their main features and comparison with other genomes. We used massively parallel 454 pyrosequencing to obtain a 2x coverage of wheat chromosome 5A. The resulting sequence assembly was used to identify TEs, genes and miRNAs, as well as to infer a virtual gene order based on the synteny with other grass genomes. Repetitive elements account for more than 75% of the genome. Gene content was estimated considering non-redundant reads showing at least one match to ESTs or proteins. The results indicate that the coding fraction represents 1.08% and 1.3% of the short and long arm respectively, projecting the number of genes of the whole chromosome to approximately 5,000. 195 candidate miRNA precursors belonging to 16 miRNA families were identified. The 5A genes were used to search for syntenic relationships between grass genomes. The short arm is closely related to Brachypodium chromosome 4, sorghum chromosome 8 and rice chromosome 12; the long arm to regions of Brachypodium chromosomes 4 and 1, sorghum chromosomes 1 and 2 and rice chromosomes 9 and 3. From these similarities it was possible to infer the virtual gene order of 392 (5AS) and 1,480 (5AL) genes of chromosome 5A, which was compared to, and found to be largely congruent with the available physical map of this chromosome.

Signaling pathways mediating the induction of apple fruitlet abscission.

Apple (Malus × domestica) represents an interesting model tree crop for studying fruit abscission. The physiological fruitlet drop occurring in this species can be easily magnified by using thinning chemicals, such as benzyladenine (BA), to obtain fruits with improved quality and marketability. Despite the economic importance of this process, the molecular determinants of apple fruitlet abscission are still unknown. In this research, BA was used to obtain fruitlet populations with different abscission potentials to be analyzed by means of a newly released 30K oligonucleotide microarray. RNAs were extracted from cortex and seed of apple fruitlets sampled over a 4-d time course, during which BA triggers fruit drop, and used for microarray hybridization. Transcriptomic profiles of persisting and abscising fruitlets were tested for statistical association with abscission potential, allowing us to identify molecular signatures strictly related to fruit destiny. A hypothetical model for apple fruitlet abscission was obtained by putting together available transcriptomic and metabolomic data. According to this model, BA treatment would establish a nutritional stress within the tree that is primarily perceived by the fruitlet cortex whose growth is blocked by resembling the ovary growth inhibition found in other species. In weaker fruits, this stress is soon visible also at the seed level, likely transduced via reactive oxygen species/sugar and hormones signaling cross talk, and followed by a block of embryogenesis and the consequent activation of the abscission zone.

Different phylogenomic approaches to resolve the evolutionary relationships among model fish species.

Comparative genomics holds the promise to magnify the information obtained from individual genome sequencing projects, revealing common features conserved across genomes and identifying lineage-specific characteristics. To implement such a comparative approach, a robust phylogenetic framework is required to accurately reconstruct evolution at the genome level. Among vertebrate taxa, teleosts represent the second best characterized group, with high-quality draft genome sequences for five model species (Danio rerio, Gasterosteus aculeatus, Oryzias latipes, Takifugu rubripes, and Tetraodon nigroviridis), and several others are in the finishing lane. However, the relationships among the acanthomorph teleost model fishes remain an unresolved taxonomic issue. Here, a genomic region spanning over 1.2 million base pairs was sequenced in the teleost fish Dicentrarchus labrax. Together with genomic data available for the above fish models, the new sequence was used to identify unique orthologous genomic regions shared across all target taxa. Different strategies were applied to produce robust multiple gene and genomic alignments spanning from 11,802 to 186,474 amino acid/nucleotide positions. Ten data sets were analyzed according to Bayesian inference, maximum likelihood, maximum parsimony, and neighbor joining methods. Extensive analyses were performed to explore the influence of several factors (e.g., alignment methodology, substitution model, data set partitions, and long-branch attraction) on the tree topology. Although a general consensus was observed for a closer relationship between G. aculeatus (Gasterosteidae) and Di. labrax (Moronidae) with the atherinomorph O. latipes (Beloniformes) sister taxon of this clade, with the tetraodontiform group Ta. rubripes and Te. nigroviridis (Tetraodontiformes) representing a more distantly related taxon among acanthomorph model fish species, conflicting results were obtained between data sets and methods, especially with respect to the choice of alignment methodology applied to noncoding parts of the genomic region under study. This may limit the use of intergenic/noncoding sequences in phylogenomics until more robust alignment algorithms are developed.

Short-term postharvest carbon dioxide treatments induce selective molecular and metabolic changes in grape berries.

Detached wine grapes ( Vitis vinifera cv. 'Trebbiano', white skinned) were treated for 3 days with 30 kPa of CO(2) and then transferred to air for an additional 9 days to partially dehydrate (about 20% weight loss). At the end of the CO(2) treatment on withering berries, total polyphenols and flavonoids were maintained in the skin, but to a more limited extent in the pulp. An induction of the proanthocyanidin synthesis appeared to be one of the responses to the treatment because both (+)-catechin and (-)-epicatechin concentrations increased in the skin. The skin and pulp of the grape berries showed different molecular responses to a high CO(2) treatment. As revealed by microarray hybridizations, 217 and 75 genes appeared differentially expressed in the skin and pulp of treated samples, respectively. Functional categorization and gene enrichment analyses pointed out that epicarp cells undergo more pronounced changes in transcript profiling at the end of the incubation period. Highly represented categories in both tissues were related to protein, stress, transcript, RNA, and hormone (ethylene, ABA) metabolism. Fermentation, CHO metabolism, and redox regulation functional categories were represented only in the skin.

PASS: a program to align short sequences.

SUMMARY: Standard DNA alignment programs are inadequate to manage the data produced by new generation DNA sequencers. To answer this problem, we developed PASS with the objective of improving execution time and sensitivity when compared with other available programs. PASS performs fast gapped and ungapped alignments of short DNA sequences onto a reference DNA, typically a genomic sequence. It is designed to handle a huge amount of reads such as those generated by Solexa, SOLiD or 454 technologies. The algorithm is based on a data structure that holds in RAM the index of the genomic positions of 'seed' words (typically 11 and 12 bases) as well as an index of the precomputed scores of short words (typically seven and eight bases) aligned against each other. After building the genomic index, the program scans every query sequence performing three steps: (1) it finds matching seed words in the genome; (2) for every match checks the precomputed alignment of the short flanking regions; (3) if passes step 2, then it performs an exact dynamic alignment of a narrow region around the match. The performance of the program is very striking both for sensitivity and speed. For instance, gap alignment is achieved hundreds of times faster than BLAST and several times faster than SOAP, especially when gaps are allowed. Furthermore, PASS has a higher sensitivity when compared with the other available programs. AVAILABILITY AND IMPLEMENTATION: Source code and binaries are freely available for download at http://pass.cribi.unipd.it, implemented in C++and supported on Linux and Windows.